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ABSTRACT 

We study the dynamical instability against bar-mode deformation of differentially 
rotating stars. We performed numerical simulation and linear perturbation analysis 
adopting polytropic equations of state with the polytropic index n = 1. It is found 
that rotating stars of a high degree of differential rotation are dynamically unstable 
even for the ratio of the kinetic energy to the gravitational potential energy of 0(0.01). 
Gravitational waves from the final nonaxisymmetric quasistationary states are calcu- 
lated in the quadrupole formula. For rotating stars of mass 1.4Mq and radius several 
10 km, gravitational waves have frequency several 100 Hz and effective amplitude 
- 5 x 10~ 22 at a distance of ~ 100 Mpc. 
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1 INTRODUCTION 



Stars ii 
ric rota 



Tohline 1987,1988; Tohline fe Hachisu 199C; Houser et al 



nature arc rotating and .subject to nonaxisymmet- 



OTirm iHsrarannss: ah ma rxmmm, w am in- 



stabilities easts only for Incompressible cicjulllbrlLim IluhFr 



199^; [Smith et al. 1995|; [Houser fe Centrella 1996|; |Picke~ 



et al. 19961 ; [Toman et al. 1998|; [mamura et al. 2000|; |New 
et al. 200C; Brown 200C|; |Liu fe Lindblom 2001]; [Liu 200 1|) 



in Newtonian gravity ( |Chandrasekhar 1969] ; Tassoul 1978). 
For these configurations, global rotational instabilities arise 
from nonradial toroidal modes e lmip (m — ±1, ±2, . . .) when 
= T/W exceeds a certain critical value. Here (p is the az- 
imuthal coordinate and T and W are the rotational kinetic 
and gravitational potential energies. In the following we will 



have shown that 0d is ~ 0.27 as long as the rotational pro- 
file is not strongly differential. In this case, once a bar has 
developed, the formation of spiral arms plays an important 
role in redistributing the angular momentum and forming a 
core-halo structure. Recently, it has been shown that 0a can 
be s maller for stars with a higher degree of different i al rota- 



focus oih the m = ±2 bar-mode since it is expected to be the 



fastest growing mode (but see Centrella et al. (zuul) with 
regard to a counter-example for extremely soft equations of 
state). 

There exist two different mechanisms and correspond- 
ing timescales for bar-mode instabilities. Uniformly rotating, 
incompressible stars in Newtonian theory are secularly un- 
stable to bar- mode deformation when > S ~ 0.14. How- 
ever, this instability can only grow in the presence of some 
dissipative mechanism, like viscosity or gravitational radia- 
tion, and the growth time is determined by the dissipative 
timescale, which is usually much longer than the dynamical 
timescale of the system. By contrast, a dynamical instabil- 
ity to bar-mode deformation sets in when > 0d — 0.27. 
This instability is independent of any dissipative mecha- 
nisms, and the growth time is determined typically by the 
hydrodynamical timescale of the system. 



et al. 1996 



tion (Centrella ct al. 2001; Tohline fc Hachisu 199C 



Liu fe Lindblom 2001; Liu 2001 ) as f3 d ~ 0.14 



Pickett 



Fo: the compressible case, determining the onset of the 



dynamical bar-mode instability, as well as the subsequent 
evolution of an unstable star, requires numerical computa- 
tions. H ydrodynamical simulations performed i n Newtonian 
theory (Tohline et al. 1985; Durisen et al. 1986; Williams & 



In such case, the formation of bars and spiral arms does not 
take place. Instead, small density perturbation is excited and 
left to be weakly nonlinear. To date, there is no report of 
dynamically unstable stars with <J 0.14. 

There are numerous evolutionary paths which may lead 
to the formation of rapidly and differentially rotating neu- 
tron stars. increases approximately as during stellar 
collapse. Also, with decreasing stellar radius, the angular ve- 
locity in the central region may be much larger than that in 
the outer region. During supernova collapse, the core con- 
tracts from ~ 1,000 km to several 10 km, and hence in- 
creases by two orders of magnitude. Thus, even rigidly rotat- 
ing progenitor stars with a moderate angular velocity may 
yield rapidly and differentially rotating neutro n stars which 
may reach the onset of dynami cal instability (Bonazzola & 
Marck 1995| ; |Rampp et al. 1998| ). Similar arguments hold for 



accretion induced collapse of white dwarfs to neutron stars 
and for the merger of binary white dwarfs to neutron stars. 
Differential rotation is eventually suppressed by viscous an- 
gular momentum transfer or magnetic braking (Baumgarte 
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2000), but within the transport timescale, the differentially 
rotating stars are subject to nonaxisymmetric instabilities. 

In this Letter, we report some of results we have re- 
cently obtained on dynamical bar-mode instabilities in dif- 
ferentially rotating stars. We pay attention to rotating stars 
with a high degree of differential rotation in Newtonian grav- 
ity using both linear perturbation analysis and nonlinear 
hydrodynamical simulation. T he detail of our study will b e 
summarized in a future paper (Shibata et al. in preparation). 
Here, we highlight a new finding. 
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2 METHOD 

To investigate the dynamical stability, we first prepare ro- 
tating stars in equilibrium. Rotating stars in equilibrium are 
modeled by the polytropic equations of state as P = Kp T 
where P, p, K and F denote the pressure, density, polytropic 
constant and adiabatic index. Here, we choose T — 2 to 
model moderately stiff equations of state for neutron stars. 
For simplicity, the profile of the angular velocity il is set as 

Q A 2 



Q = 



+ A 2 ' 



(1) 



where A is a constant, Qo the angular velocity at the sym- 
metric axis, and m cylindrical radius (a/ 1 x 2 + y 2 ). A controls 
the steepness of the rotational profile: For smaller A, the pro- 
file is steeper and for A — » oo, the rigid rotation is recovered. 
In the limit of A going to 0, this rotation law gives uniform 
specific angular momentum, which is the limiting law al- 
lowed by axisymmetric stability analyses. (Here, however, 
A = is prohibited in this situation because angular veloc- 
ity diverges along the rotation axis.) This is also the rotation 
law that has been often used in studies of nonaxisymmetric 



instabilities in tori and an nuli ( Papaloizou & Pringle 1984 
Goodman fc Narayan 1988| ; |Andalib et al. 1997[ ). In this Let- 



ter, we report the results for A = A/R = 0.3 and 1 where R 
is the equatorial radius of rotating stars. 

In a linear perturbation analysis, we derive the linear 
perturbation equations of the hydrodynamic equations and 
substitute the functions in the form f(r, 0) e lm ¥> _la '* f or p er _ 
turbed hydrodynamic variables. As a result, the problem re- 
duces to the eigen value problem for determining a complex 
value of u and the corresponding eigen functions of the per- 
turbed quantities. The details of the numerical method w ill 
be reported in future (Karino & Eriguchi in preparation). 



In the hydrodynamical simulation, we initially superim- 
pose a density perturbation of the bar-mode in the form 

,,2 



Sp- 



Po- 



y 



R 2 



(2) 



where po denotes the axisymmetric configuration and S con- 
stant. Namely, we focus on a fundamental mode (f-mode) for 
which there is no node. In this work, we choose 5 = 0.1, and 
the velocity is left to be unperturbed at t = 0. The growth 
of a bar-mode can be followed by monitoring the distortion 
parameter 



rj ; 



+ Iy 



where hj denotes the quadrupole moment 
lij = / d 3 xpx'x J ] . 



(3) 



(4) 



Figure 1. \r)\ as a function of for A = 1 with C a = 0.255 
(solid line) and 0.305 (dashed line). The dotted line denotes the 
growth rate of the unstable mode. 



Simulations were p erformed using a 3D numerical code in 
Newtonian gravity ( |Shibata et al. 1997 ). We adopt a fixed 
uniform grid with size 141 x 141 x 141 in x — y — z, which 
covers an equatorial radius by 50 grid points initially. To 
confirm that the results depend weakly on the resolution, 
we also perform test simulations with size 71 x 71 x 71 (i.e., 
the grid spacing becomes twice larger) for several selected 
cases. We assume a reflection symmetry about the equatorial 
plane. Since many rotating stars we pick up have flattened 
configuration, we set the grid spacing of z half of that of x 
and y. 



3 NUMERICAL RESULTS 

First we show results for rotating stars of moderately large 
differential rotation with A = 1. In Fig. 1, we show |?7| as a 
function of time for C a = R P /R = 0.255 and 0.305, where 
R p is the polar radius. It is found that for C a = 0.255, 
the rotating stars are dynamically unstable, while for C a = 
0.305, they are stable. In Fig. 2, we also display w r /f2o and 
u)i/Q.a as a function of /3. Here, u) r was derived by carrying 
out the Fourier transformation of r)(t), and Ui was measured 
using the growth rate of the peaks of \r)\ (see dotted line in 
Fig. 1). We also compute these angular frequencies by a 
linear perturbation analysis (solid lines). We find that two 
results independently obtained agree fairly well. The linear 
analysis also indicates that this unstable mode is likely to 
be the f-mode (of no node in the eigen function). 

We find that at C a ~ 0.3, the dynamical stability 
changes. In terms of 13, the stability changes at f3 ~ 0.26 
for A = 1. Thus, (3d is not extremely smaller than that fo r 
rigidly rotating incompressible stars ( Chandrasekhar 196£ ) 



in the case when the degree of the differential rotation is not 
very high. 

The situation drastically changes for A — 0.3. In Fig. 
3, we show |?7| as a function of time for C a = 0.405, 0.605, 
and 0.805 obtained by numerical simulations. Note that for 
C a = 0.805, P » 0.04. We find that in every model, non- 
axisymmetric dynamical instabilities take place. In Fig. 4, 
we also show u} r /Q.o and Wi/fio a s a function of j3. Since 
several growing modes seem to be excited simultaneously 
for C a = 0.605 and 0.805, we plot upper and lower bounds 
for Ui. Results by numerical simulations agree again with 
those by linear perturbation analysis. The linear analysis 
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Figure 2. uj r /flo and Wi/ilo as a function of /3 for A = 1. The 
solid curves and solid circles denote the results by linear pertur- 
bation analysis and numerical simulation, respectively. We choose 
results for C a = 0.205, 0.255, 0.305 and 0.355. 



C =0.405 




Figure 3. |?j| as a function of Qyt for A = 0.3 with C a = 0.405 
(upper) 0.605 (middle), and 0.805 (lower). The dotted lines de- 
note the growth rate of the unstable modes. For C a = 0.605 and 
0.805, we plot two dotted lines which indicate the upper and lower 
bounds of the growth rates. 



also indicates that even for j3 fa 0.03, the rotating stars are 
dynamically unstable. 

In Fig. 5, we display the contour line and velocity vec- 
tors for C a — 0.605 at selected time slices as an example. 
We find that the dynamical instability does not excite large 
spiral arms and fragmentations. Rather, it saturates in a 
weakly nonlinear stage at which |7j| ~ 0.2, to form a quasis- 
tationary ellipsoid. We note that this feature is independent 
of C a as long as A = 0.3. 

Since the rotating stars of C a Ss 0.85 change to non- 
axisymmetric objects, they can be sources of gravitational 
waves. In Fig. 6, we show the gravitational waveforms along 
the z— axis (h+ tX ) and energy luminosity (E) as a function 
of time for C a = 0.605. We note that the results for other 
models are qualitatively the same as that in Fig. 6. Here, 




Figure 4. The same as Fig. 2, but for A = 0.3. For C a > 0.605, 
upper and lower bounds of aj; for results in numerical simulations 
are shown. 



gravitational waves are calculated in the quadrupole formula 
( |Misner et al. 1975j ), and we define 



■21 



r 



and E : 



5 



(3)^(3) 



(5) 



(3) 



where iy = d 2 Iij/dt 2 , 4-y = hj — SijI kk /3, 4- 
d 3 4rij I 'dt 3 , and r is the distance from a source to a detec- 
tor. Here, we use the units c = G = 1. Since the outcome 
of the dynamical instability is a quasistationary ellipsoid, 
the amplitude and luminosity of gravitational waves settle 
down to a nearly constant value as rh+, x ~ 0.2(M 2 / R) and 
E ~ 0.01(M/i?) 5 . 



4 DISCUSSIONS 

Using the results shown in Fig. 6, we can estimate the ex- 
pected amplitude of gravitational waves from the nonax- 
isymmetric outcome formed after the dynamical instability 
saturates. Here, we pay particular attention to protoneu- 
tron stars likely formed soon after the supernovae of mass 
~ IAMq and radius several 10 km. For the maximum am- 
plitude of gravitational waves h ~ 0.2M 2 /Rr, and the en- 
ergy luminosity E ~ 0.01 (M/R) , the emission timescale of 
gravitational waves can be estimated as 



100/3' ( ^ 
E \M 



M. 



(6) 



where we set T = f3'M 2 /R and 0.04 & f3' <, 0.4 for A = 
0.3 and 0.3 < C a < 0.8. The characteristic frequency of 
gravitational waves is denoted as 



2tt 



800Hz 



/ 15M \ 
V R J 



3/2 



M 



1.4M, 



(7) 



where ^ = u r ^jR 3 /M ~ 2 (see Figs. 4 and 7). Since the 
cycles of gravitational wave-train N are estimated as N = 
fr, the effective amplitude of gravitational waves h e g = 
y/lfh is 



h e 



5 x 10" 



1/2 



£_ 

o.i 



1/2 



15M J 



1/4 
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Figure 6. Gravitational waves for A = 0.3 and C a = 0.605. 
Upper figure: h-y. (solid curve) and h x (dotted curve) in units of 
M 2 /R as a function of Qq*- Lower figure: E in units of (M/R) 5 . 




X / R 



Figure 7. Uq/t/m/R? as a function of /3. The solid circles de- 
note the dynamically unstable stars we found. 



Figure 5. Snapshots of density contour lines and velocity vectors 
in the equatorial plane at t = (left) and Qot PS 390 (right) 
for A = 0.3 and C a = 0.605. The contour lines are drawn for 
p/pmax = j/10 (J = 1 ~ 9) and 0.95 (solid curves) and for 
p/Pmax = 0.01 (dotted curve), where p m ax denotes the maximum 
density. 



M \( lOOMpc \ 
1-4M J \ r ) 



(8) 



([Thorne 1987| ; |Lai fc Shapiro 199E| ; [Liu fc Lindblom 2001 



Liu 2001). Thus, gravitational waves from protoneutron 



stars of a high degree of differential rotation and of mass 
~ IAMq and radius several 10 km in the distance of ~ 100 
Mpc can be a source for laser interferometric detectors such 
as LIGO ( Thorne 1995 ). We emphasize that / and h a a de- 
pend weakly on f3 (or j3') for A = 0.3. This implies that even 
if the star is not rapidly rotating, the dynamical instability 
can set in and as a result, differentially rotating stars can 
emit gravitational waves of a large amplitude. 

To summarize, we have studied dynamical bar-mode in- 
stability of differentially rotating stars focusing on the f- 
mode. We have found that rotating stars of a high degree of 
differential rotation are dynamically unstable against non- 
axisymmetric deformation even for (3 <C 0.27. 



It is worthy to note that the real parts of the eigen fre- 
quencies do not vanish but approach to finite values as the 
value of /3 decreases, i.e. in the spherical limit. This behav- 
ior is totally different from those of the r-modes (see e.g. 
Karino et al. 2001), several self-gravity induced instability 
modes of slender tori or annu l i such as I-modes an d J-modes 
( |Goodman fc Narayan 1988 ; Andalib et al. 1997), or shear 



instability modes (or P-modcs) such as Papaloizou -Pringle 
instability for toroids ( 3 apaloizou fc Pringle 198^ ) and for 
spheroids (Luyten 199C). Therefore, we have identified the 
unstable modes we find in this paper as the f-mode. 

The physical mechanism for the onset of dynamical in- 
stabilities found in this Letter may be explained in the fol- 
lowing manner: For a small value of A, non-spherical de- 
formation of a stellar configuration is significant around the 
rotational axis although the distant part from the rotational 
axis is almost spherical. Thus, with decrease of the value of 
C a , most rotational energies are confined to the region near 
the rotational axis. Increase of the rotational energy can be 
much larger than those of the gravitational and internal en- 
ergies, because the overall shape cannot be much different 
from that of a sphere. As a result, the total energy becomes 
large with small values of C a - However, if the rotational en- 
ergy exceeds a certain amount in the region near the rota- 
tional axis, there may exist other equilibrium configurations 
with lower total energies, as in the case of bifurcation of 
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t he Jacobi ellipsoidal sequence from the Maclaurin sequence 
( Chandrasekhar 1969). In fact, the total energy of the Jacobi 



ellipsoid is lower than the Maclaurin spheroid of the same 
mass and the angular momentum, if the rotational energy 
exceeds a certain criterion. 

Once the instability sets in, the axisymmetric star be- 
gins to change its shape to a nonaxisymmetric configu- 
ration. In such case, the angular velocity in the region 
near the rotational axis, ~ Qo, is very large (larger than 
y/ M/R 3 ). It implies that, despite of the large rotational 
energy, the angular momentum cannot be very large be- 
cause T ~ f(x 2 + y 2 ) x f^dm and J ~ f(x 2 + y 2 ) x Q dm, 
where J and dm are the angular momentum and the mass 
element, respectively. Therefore, even if the nonaxisymmet- 
ric mode grows due to the dynamical instability, the con- 
figuration cannot be highly elongated or form spiral-arm 
structures; i.e., there exists a saturation state of a small 
nonaxisymmetric deformation. For the bar-mode instability 
of the Maclaurin spheroids, the rotational energy is "con- 
fined" to the outer part of the configuration by making the 
shape very flat and so the angular momentum can be also 
large, even th ough the angular velocity is not extremely large 
(< y/M/R 3 ). Therefore, once a bar-mode dynamical insta- 
bility sets in, configurations can be considerably different 
from the original spheroidal shapes. 

It is possible that other non-fundamental modes have 
larger u)t than the f-mode has. However, the present study 
shows that there exists at least one mode which induces 
the nonaxisymmetric deformation. Such rotating stars sub- 
sequently form nonaxisymmetric structures and hence can 
be a source of laser interferometric gravitational wave de- 
tectors. 
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